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We test the rehabiUty of the the Polyakov-loop extended Nambu-Jona-Lasinio (PNJL) model, comparing the 
model result with the lattice data at nonzero imaginary chemical potential. The PNJL model with the vector-type 
four-quark and scalar-type eight-quark interactions reproduces the lattice data on the pseudocritical temperatures 
of the deconfinement and chiral phase transitions. The QCD phase diagram in the real chemical potential region 
is predicted by the PNJL model. The critical endpoint survives, even if the vector-type four-quark interaction is 
taken into account. 

PACS numbers: 11.30.Rd, 12.40.-y 



I. INTRODUCTION 

Quantum Chromodynamics (QCD) is a remarkable theory. 
It is renormalizable and essentially parameter free. QCD ac- 
counts for the rich phenomenology of hadronic and nuclear 
physics. Thermodynamics of QCD is also well defined. Nev- 
ertheless, it is not well known because of its nonperturbative 
nature. In particular, QCD phase diagram is essential for un- 
derstanding not only natural phenomena such as compact stars 
and the early universe but also laboratory experiments such as 
relativistic heavy-ion collisions. 

Unfortunately, quantitative calculations of the phase dia- 
gram from first-principle lattice QCD (LQCD) have the well 
known sign problem when the chemical potential (/i) is real; 
for example, see Ref. 1 1 ] and references therein. So far, sev- 
eral approaches have been proposed to circumvent the diffi- 
culty; for example, the reweighting method |2], the Taylor 
expansion method 13|1 and the analytic continuation to real 
chemicalpotential (/zr) from imaginary chemical potential 
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fection. 



However, those are still far from per- 



As an approach complementary to first-principle lattice 
QCD, we can consider effective models such as the Nambu- 
Jona-Lasinio (NJL) model |i,[l3,[ni[ll[ll[ll[Il[ll 
[nil and the Polyakov-loop extended Nambu-Jona-Lasinio 
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35[l . The NJL model describes the chiral 



symmetry breaking, but not the confinement mechanism. The 
PNJL model is designed fid] to make it possible to treat the 
Polyakov loop as well as the chiral symmetry breaking. 

In the NJL-type models, the input parameters are deter- 
mined at fi = and T > 0, where T is temperature. It is 
then highly nontrivial whether the models predict properly dy- 
namics of QCD at finite /iR. This should be tested from QCD. 
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Fortunately, this is possible in the fii region, since lattice QCD 
has no sign problem there. The canonical partition function 
Zc (n) with real quark number n is the Fourier transform of 
the grand-canonical one Zqc{9) with 9 = Hi/T l[36ll : 
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Thus, the thermodynamic potential of QCD, f^QcriiS) — 
— T1ii(Zgc(^)), at finite 9 includes all dynamics at real n 
and hence at finite /i^. Therefore, the reliability of effective 
models at finite fi^ can be tested in the fii region. 

Roberge and Weiss found ll36ll that QCD has a periodicity 
^QCDiO) = X?qcd(6' + 27rfc/3), showing that I2qcd(6' + 
27rfc/3) is transformed into ^2qcd(^) by the Z3 transforma- 
tion with integer k. This means that QCD is invariant under a 
combination of the Z3 transformation and a parameter trans- 
formation 9^9 + 2kn/3 Hill, 

q Uq, A^^UAM'^ -i/g{dM)U-\ 

9 + 27rfc/3, (2) 

where U{x, t) are elements of SU(3) with U{x, (3 = l/T) = 
exp(— 2i7rfc/3)[/(a;, 0) and q is the quark field. We call this 
combination the extended Z3 transformation. Thus, i7QCD(^?) 
has the extended Z3 symmetry, and hence quantities invariant 
under the extended Z3 transformation have the RW periodic- 
ity ll34i Issil . At the present stage, the PNJL model is only a 
realistic effective model that possesses both the extended Z3 
symmetry and chiral symmetry jHHH]. This property makes 
it possible to compare PNJL with lattice QCD quantitatively 
in the /ij region. If the PNJL model succeeds in reproducing 
the lattice data, we may think that the PNJL model will pre- 
dict, with high reliability, the QCD phase structure in the /iR 
region. 

The extended Z3 symmetry in QCD is a remnant of the Z3 
symmetry, namely the confinement mechanism, in the pure 
gauge system. The extended Z3 symmetry appears as the RW 
periodicity in the /ij region and implicitly affects dynamics in 
the /iR region. Actually, the mechanism largely shifts the crit- 
ical endpoint [ 10] toward higher T and lower fi than the NJL 
model predicts 1 21, 25, 31]. In contrast, the vector-type four- 
quark interaction Gv{qjfj,q)'^ largely moves the critical end- 
point in the opposite direction UA. Im I21I Isill . if it is newly 
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added to the NJL and PNJL models. Thus, it is essential to de- 
termine the strength of the coupling of the vector-type in- 
teraction, although the interaction is often ignored in the NJL 
and PNJL calculations. 

In the relativistic meson-nucleon theory iIjtIi . the repulsive 
force mediated by vector mesons is essential to account for the 
saturation property of nuclear matter. Using the auxiliary field 
method, one can convert quark-quark interactions to meson- 
quark interactions; for example, see Refs. ITtI [ssi [39ll and 
references therein. In the hadron phase, quarks have a large ef- 
fective mass as a result of spontaneous chiral symmetry break- 
ing, and then nucleons can be considered to be formed by such 
three heavy quarks, i.e. three constituent quarks. It is then nat- 
ural to think that there exists the correspondence between the 
meson-nucleon interactions and the quark-quark interactions. 
In this sense, it is very likely that the vector-type four-quark 
interaction is not negligible and even significant in particular 
at a finite quark-density region corresponding to the nuclear 
saturation density. In the previous work |35], we have pro- 
posed that the strength of Gv can be determined from lattice 
data on the chiral phase transition in the /ii region. 

In this paper, we consider two-flavor QCD and show the 
reliability of the PNJL model, quantitatively comparing the 
model result with lattice data in the jii region. The model pa- 
rameters except Gv are fixed by the measured pion mass and 
decay constant at /i = T = and lattice data [4^ I4lil42ll at 
T > and /i = 0. The PNJL calculation with no vector-type 
interaction well reproduces lattice data iS |^ on the pseudo- 
critical temperature Tc(^) of the deconfinement phase transi- 
tion, but not on the pseudocritical temperature Tc{(j) of the 
chiral phase transition near 9 = vr/S. The strength of Gv is 
fitted so as to reproduce the latter data. The primary result 
of the lattice simulations is that Tc{(l>) coincides with rc(cr), 
within numerical errors, in the entire region of 6* [sj] . The 
PNJL model with the vector-type interaction can reproduce 
this property. Finally, we quantitatively predict the phase di- 
agram in the /iR region by using the PNJL model with the 
parameter set justified in the /ij region. This sort of model 
predictions are quite important before doing heavy lattice cal- 
culations with large lattice size in the /ii region. 

In section II, the PNJL model is explained simply. In sec- 
tion III, we test the PNJL model in the /ij region and determine 
the strength of Gv. Finally, we predict the phase diagram in 
the /xr region. Section IV is devoted to summary. 



interaction. Later, we will add the vector-type four-quark in- 
teraction ITc l, [T^ [l6i I35I1 and the scalar-type eight-quark in- 
teraction 1I15I [m I34I] to the PNJL Lagrangian. The Polyakov 
potential U, defined in dHJ, is a function of the Polyakov loop 
and its Hermitian conjugate <!>*, 



with 
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L{x) =V exp 



drA. 
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where V is the path ordering and A4 — iAo. In the chiral limit 
(mo ~ 0), the Lagrangian density has the exact SU{Nf)i, x 
SU{Nf )R X C/(l)v X SU{3)c symmetry. 

The temporal component of the gauge field is diagonal in 
the flavor space, because the color and the flavor space are 
completely separated out in the present case. In the Polyakov 
gauge, L can be written in a diagonal form in the color 
space I20II : 

^ ^ p,:/3(03A3+08A8) ^ diag(e*'^'^% 6*'^'^% e''^*'^), (6) 

where (pa = (ps + (ps/VS, (pt = ~(p3 + (ps/VS and (pc = 
— {(pa + (pb) = —2(ps/VS- The Polyakov loop ^ is an exact 
order parameter of the spontaneous Z3 symmetry breaking in 
the pure gauge theory. Although the Z3 symmetry is not an 
exact one in the system with dynamical quarks, it still seems 
to be a good indicator of the deconfinement phase transition. 
Therefore, we use (p to define the deconfinement phase transi- 
tion. 

Making the mean field approximation and performing the 
path integral over quark field, one can obtain the thermody- 
namic potential i7 (per volume). 



f2 



2Nf 



d^P r 



(27r)3 



- In [1 + 3{<P + <?*e-/5^^-(p))e-'5^"(P) + e-3/5S-(p)] 
P 

- In [1 + 3{<P* + <Z.e-'5^^+(p))e-'^^^(P) + g-^/^^^^P)]' 

P 



(7) 



II. PNJL MODEL 

The two-flavor PNJL Lagrangian is 

£ ^q{ij^D'' - mo)q 

+ Gs[{qqf + {qtj5fqf]~U{'P[Al<P[A]*,T), (3) 

where q denotes the two-flavor quark field, toq does the cur- 
rent quark mass, and I?" = d'' + iA'^ — ifiSg. The field A^ is 
defined as A" = 6^gA°^ with flie gauge field A^, the Gell- 
Mann matrix Aa and the gauge coupling g. In the NJL sector, 
Gs denotes the coupling constant of the scalar-type four-quark 



where, <t = (qq), Sg = — SGgCr, M = mo + Ss, Um — GgCr^, 
E{p) = mdE^ip) = E{p)±n = E{p)±i9/f3. 

In (|7]i, only the first term of the right-hand side diverges. It is 
then regularized by the three-dimensional momentum cutoff 
A ^ g4J. We use U of Ref. [25] that is fitted to a lattice 
QCD simulation in the pure gauge theory at finite T ll43l 01 '■ 



U = T'' 



a{T) 



b{T) ln(l - Q^P 



*3\ 



3 )^ 



, (8) 



a{T) = ao + ai(^) + a2(f b{T) ^ 63 (y)' (9) 
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where parameters are summarized in Table I. The Polyakov 
potential yields a first-order deconfinement phase transition at 
T = To in the pure gauge theory. The original value of Tq 
is 270 MeV evaluated by the pure gauge lattice QCD calcula- 
tion. However, the PNJL model with this value of Tq yields 
somewhat larger value of the transition temperature at zero 
chemical potential than the full LQCD simulation ll40ll4lll42ll 
predicts. Therefore, we rescale Tq to 212 MeV; the detail will 
be shown in subsection llll Al 



ao 


ai 


02 


b3 


3.51 


-2.47 


15.2 


-1.75 



TABLE I: Summary of the parameter set in the Polyakov sector used 
in Ref. flB]. All parameters are dimensionless. 



The explicit 9 dependence appears only through the factor 
e^*^ in (fT4l i. Hence, the stationary conditions (fTOb show that 
X — X{e^'^^). Inserting the solutions back to (fT4l) . one can 
see that H = Q{e^'^). Thus, X and Q have the RW periodic- 
ity, 

X{e+=^)=X{0), and f2(9+^) = f2ie), (15) 

o o 

while the Polyakov loop ^ and its Hermitian conjugate 
have the properties 

^0+—)* = e*2'^'=/3<?(6')*. (16) 



The variables X ^ (l>, <P* and a satisfy the stationary con- 
ditions. 



^^2/^X = 0. 



(10) 



The solutions of the stationary conditions do not give the 
global minimum f2 necessarily. There is a possibility that they 
yield a local minimum or even a maximum. We then have 
checked that the solutions yield the global minimum when the 
solutions X{9) are inserted into (|7]i. 

The thermodynamic potential J? of Eq. (jT) is not invariant 
under the Z3 transformation. 



^9) ^>(6')e-*2^'=/3 , <P{9y <P{9) 



(11) 



although of ([8]) is invariant. Instead of the Z3 symmetry, 
however, f2 is invariant under the extended Z3 transformation. 



<P{9y 



${9ye' 



<P{9) <P{9)e- 
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ni. NUMERICAL RESULTS 
A. Thermal system with no chemical potential 

First, we consider the thermal system with no chemical 
potential to determine the parameters, mo, Gs, A and Th of 
the PNJL model. In the lattice calculations H liH l42ll . 
the pseudocritical temperature Tc{(j) of the crossover chiral 
phase transition coincides with that Tc{'P) of the crossover 
deconfinement one within 10 % error: Tc(o') w Tc{'P) ~ 
173±8MeV||4lll. 

The parameter set, A = 631.5 MeV, Gs = 5.498 [GeV-^] 
and mo = 5.5 MeV, can reproduce the pion decay constant 
= 93.3 MeV and the pion mass = 138 MeV at 
T = fi = lfl6ll . and keeps a good reproduction also at fi- 
nite T ll25ll . We then adopt these values for A, Gs and mo. We 
adjust To so that the PNJL calculation can reproduce the lat- 
tice result Tc(<P) = 173 MeV; the value is To =212 MeV. The 
parameter set thus determined is shown as set A in Table HIl 



This is easily understood as follows. It is convenient to intro- 
duce the modified Polyakov loop 'F = e'^^ and 'F* = e^'^^* 
invariant under the transformation ( fT2] l. The extended Z3 
transformation is then rewritten into 



W{9)* 

and i7 is also into 



-> •f(6l)*, 



(13) 



i7 : 



2Nf 



(2^ 



3E{p) 



I In [1 + 3tf'e-^^(P) 



f3 



+ i In [1 + 3.Z/*e-'3^(P) + 3<Fe-^mp)e-0f^B 



g-3/3£;(p)g-/3AiE 



(14) 



where /3/iB — 3/3 /i = 3i9. Obviously, J? is invariant under the 
extended Z3 transformation (fT3] i. since it is a function of only 
extended Z3 invariant quantities, e'^*^ and X{— ^,^*,a). 



set 




Gs8 


Gv 


A 


5.498GeV~^ 








B 


4.673GeV-2 


452.12GeV-* 





C 


4.673GcV"^ 


452.12GeV"* 


4.673GeV-2 



TABLE II: Summary of the parameter sets in the PNJL calculations. 
The parameters A, mo and To are common among the three sets; 
A = 631.5 MeV, mo = 5.5 MeV and To = 212 MeV. 



Figure [T] shows the chiral condensate a normalized by 
Co = o-\T=o,fj,=o and the absolute value of the Polyakov loop 
<P as a function of T/T^. In this paper Tc is always taken to 
be 173 MeV. The green curves represent the PNJL results of 
parameter set A, where do = —0.0302 [GeV''] in this case. 
Lattice QCD data [40, 41, 42] are also plotted by cross sym- 
bols with 10 % error bar; a and measured as a function 
of r/Tc in Refs. ® liHlil have only small eiTors, but we 
have added 10 % error that the lattice calculation |41] has in 
determining Tc. For |<?| the PNJL result (green solid curve) 
reasonably agrees with the lattice one (x). For a, however. 
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the PNJL result (green dashed curve) considerably overshoots 
the lattice data (+). 
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Fig. 1: Chiral condensate a normalized by aiT = 0, ^ = 0) and 
the absolute value of the Polyakov loop <P. The blue (green) curve 
represents the PNJL result of parameter set B (A) with (without) the 
scalar-type eight-quark interaction; a (|^|) is denoted by the dashed 
(solid) curve. Lattice data (+) on a are taken from Ref. |40] and 
those (x) on \$\ are from Ref. |42]. The lattice data are plotted 
with 10 % error bar, since lattice calculations have 10 % error in 
determining i4lll . 

Figure |2] represents results of the PNJL calculations for 
chiral and Polyakov-loop susceptibilities, Xa and x-i lUill- 
Peak positions of Xa and xs> show Tj^a) and Tc(^), respec- 
tively. The PNJL results (green curves) of parameter set A 
give Tc{a)/Tc = 1.25 and Tc(<?)/Tc = 1, while the lat- 
tice simulations yield Tc{a)/Tc = 1 ± 0.05 and Tci^)/Tc = 
1 ± 0.05. The PNJL results are consistent with the lattice ones 
for rc(<P), but not for Tc(cr). 
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Fig. 2: T dependence of chiral and Polyakov-loop susceptibilities, 
Xct (right scale) and x-p (left scale). The blue (green) curve represents 
the PNJL result of parameter set B (A) with (without) the scalar-type 
eight-quark interaction; Xc (X*) denoted by the dashed (solid) 
curve. The region between two vertical gray lines T = (1 ± 0.05)Tc 
is the prediction of lattice calculations I41i1 . 

Now we introduce the scalar-type eight-quark interac- 
tion I1I61I . 

Gssiiqqf + iqi-f^Tqff, (17) 



since the difference Tc{a) — Tc{'P) is reduced by the interac- 
tion H. 

Since and Af^ calculated with PNJL depend on the 
strength of Gss, for each value of Gss the strength of Gs 
is re-adjusted so as to reproduce the measured values = 
93.3 MeV and = 138 MeV. As Gss increases from zero, 
Tc(cr) calculated with PNJL decreases toward Tc = 173 MeV. 
When Gs8 = 452.12 GeV^^ the ratio T^ia)/T^ becomes 
1.05 and hence consistent with the corresponding lattice re- 
sult within 10 % error We adopt this strength. This parameter 
set is shown as set B in Table HIl As shown in Fig.[T] the PNJL 
results (blue curves) of parameter set B well reproduce the 
lattice results for both the chiral condensate and the Polyakov 
loop. 



B. Thermal system with imaginary chemical potential 

In this subsection, we consider the thermal system with fi- 
nite imaginary chemical potential and compare the PNJL re- 
sult with the lattice data fldSl) in which the lattice size is 
8'' X 4 and the two-flavor KS(Wilson) fermion is considered. 

First, we analyze the deconfinement phase transition. Since 
the eight-quark interaction hardly changes the Polyakov loop, 
we do the PNJL calculation with parameter set A. Figure [3] 
presents T dependence of the Polyakov-loop susceptibility X3> 
in three cases of 6* = 0, 0.56 and 0.96; each case is distin- 
guished by using different colors. For each 6, the PNJL re- 
sult (solid curve) reproduces the corresponding lattice result 
(crosses) in its peak position. Thus, the PNJL results are con- 
sistent with the lattice ones for the pseudocritical temperature 
of the crossover deconfinement phase transition. 




0.9 1 1.1 1.2 



T/Tc 

Fig. 3: T dependence of the Polyakov-loop susceptibilities in three 
cases of ^ = 0, 0.56 and 0.96; each case is distinguished by using 
different colors. The solid curves represent the PNJL results of set 
A (right scale). Lattice data shown by crosses (left scale) are taken 
from Ref. [M- 

Figure |4] presents the phase diagram of the deconfinement 
phase transition in the 6-T plane, where 6 is divided by 7r/3 
and T is normalized by Tc = 173 MeV. Lattice data [8] mea- 
sured as a function of T/Tc have only small errors, as shown 
by thick error bars in Fig. |4] This is an error bar in the case 
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that lattice calculations have no error in Tc- However, the lat- 
tice calculation [41 J has about 10 % error in determining T^., 
as mentioned in subsection IIIIAI This 10 % error should be 
added to the original small error; this 10 % error will be shown 
later in Fig. |6] The PNJL result (solid curve) of set A agrees 
with the lattice one (crosses) within the error bars. The phase 
diagram has a periodicity of 27r/3 in 9. This is called the 
Roberge and Weiss (RW) periodicity ||36ll . The phase diagram 
is also 9 even, because so is xs>- On the dot-dashed line going 
up from an endpoint (6'r{Wi Trw) — (tt/S, l.GOTc), the quark 
number density n and the phase (j) of the Polyakov loop are 
discontinuous in the PNJL calculations ll34ll35fl . This is called 
the RW phase transition line. The lattice data 10, Si on 
also discontinuous on the line, as shown later in Fi 
the PNJL result is consistent with the lattice results 
for the location of the RW phase transition line. 




H 
H 




Fig. 4: Phase diagram on the 9-T plane. The solid curve represents 
the deconfinement phase transition, while the dot-dashed line does 
the RW phase transition predicted by the PNJL calculation with set 
A. Lattice data are taken from Ref. ISl. 



The lattice simulations iJ, |8|] point out that Tc(cr) agrees 
with T'c(^) within numerical errors in the entire region < 
9 < 27r/3. We then take the case of 9 = 7r/3 to consider this 
point. It is predicted by the lattice simulations that Tc((t) and 
Tci'P) are located in the region between two vertical gray lines 
of Fig- E] Panel (a) shows cr and as a function T/Tc and 
panel (b) does Xa and as a function T/Tc- The green (blue) 
curves represent results of the PNJL calculations with set A 
(B). The eight-quark interaction hardly shifts the peak position 
of xs>, i-C- Tci'P), from the value l.GOTc. The peak position is 
consistent with the lattice result shown by the region between 
two vertical gray lines. In contrast, the eight-quark interaction 
largely shifts the peak position of Xa, i-C- 7c((t), from 1.53Tc 
to 1.24Tc, but the shifted value still deviates from Tc{<P) = 
(1.1 ± 0.05)rc, that is, the lattice data near 6* = 7r/3 1] 
shown by the region between two vertical gray lines. 

In order to solve this problem, we introduce the vector-type 
four-quark interaction 



(18) 



structure in the real chemical potential region is quite sensi- 
tive to the strength of the coupling Gv. It is then important to 
determine the strength, but it has not been done yet. Since the 
vector-type interaction does not change the pion mass and the 
pion decay constant at T = /i = and the chiral condensate 
and the Polyakov loop at T > and /i 0, we can simply 
add the interaction to set B. As Gv increases from zero, Tc{<j) 
goes down toward Tc{<P), while Tc(^) moves little. When 
Gv = 4.673 GeV~^, Tc(cr) gets into the region between the 
vertical gray lines. We adopt this strength of Gv This set is 
shown as set C in Table HIl 
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D 



and add it to the PNJL Lagrangian C; see Ref. OSp for the de- 
tail of this formulation. As mentioned in Ref. IBSll . the phase 



T/Tc 



Fig. 5: T dependence of (a) the normalized chiral condensate and 
the absolute value of the Polyakov loop and (b) the susceptibilities 
Xcr (right scale) and (left scale) at 6 = tt/S. In panel (a), a 
(|<P|) is denoted by the dashed (solid) curve. In panel (b), Xo- (x*) 
denoted by the dashed (solid) curve. The PNJL calculations are done 
with three parameter sets of A, B and C and these are distinguished 
by using different colors, green, blue and red, respectively. 

Figure |6] shows the phase diagram of the chiral phase tran- 
sition determined by Tc{a). Green, blue and red curves are 
results of the PNJL calculations with sets A, B and C, respec- 
tively. In the entire region < 9 < 2tt/3, the eight-quark in- 
teraction moves Tc((t) down from the green dashed curve (set 
A) to the blue one (set B). However, the blue dashed curve 
still overshoots the lattice result (symbols) with 10 % error 
near 9 = 7r/3. The vector-type interaction makes the blue 
dashed curve go down to the red one (set C) that is consistent 
with the lattice result |8]. Thus, the PNJL calculations with 
set C can reproduce the lattice result 10, [Ml that Tc (cr) coin- 
cides with Tc (<?) within numerical errors in the entire region 
< 61 < 27r/3. 
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Fig. 6: Phase diagrams of the chiral phase transition in the imagi- 
nary chemical potential region calculated with three parameter sets 
are presented by dashed curves; green, blue and red ones are results 
of the PNJL calculations with set A, B and C, respectively. Lattice 
data are shown with 10% error that Tc has |l4l|]. The decon- 
finement phase transition curves (solid curves) are also shown for 
comparison. 



Figure |7|a) shows 6 dependence of the phase (p of for 
four cases of T/Tc = 0.97, 1.01, 1.04 and 1.10; each case 
is distinguished by using different colors. The PNJL results 
(curves) well simulate the lattice data |4, 8] (symbols). It is 
found from both the results that (j> is continuous at = 7r/3 in 
the low-T side T < Trw = 1.09Tc, but it is discontinuous at 
6* = 7r/3 in the high-T side T > Trw Hence, the RW phase 
transition takes place at T > Trw = l.OOTc and 6 = tt/3. 

Figure |7tb) shows T dependence of (j> for five cases of 
6»/(7r/3) = 0, 0.4, 0.8, 1.0 and 1.2. The PNJL results 
(curves) well reproduce the lattice data 13- Hi (symbols). For 
9 < tt/3 the phase <j) tends to zero as T increases, while for 
6 > 7r/3 it does to — 27r/3 as T increases. When 9 — tt/3, 
the RW phase transition takes place at T > Trw = 1.09Tc 
and then the phase (j> is singular there, so that the pink line 
terminates at T = Trw- In the high-T limit, the region (I) 
-7r/3 < 9 < n/3 has </) = and the region (II) tt/3 < 9 < tt 
does (p — — 27r/3. Thus, the region (II) is a Z3 image of the 
region (I), and the region (III) tt < 61 < 57r/3 is another Z3 
image of the region (I). 



C. Thermal system with real chemical potential 

In this subsection, we predict the phase diagram in the real 
/i region by using the PNJL model. In Fig. [8] panels (a)-(c) 
represent results of the PNJL calculations with sets A, B and 
C, respectively. Panel (c) is the most reliable result, since the 
PNJL result of set C is consistent with the lattice result 
in the imaginary chemical potential region. Comparing the 
three panels, we find that the vector-type four-quark interac- 
tion and the scalar-type eight-quark interaction give sizable 
effects on the phase structure. In particular for the critical 
endpoint E, the eight-quark interaction shifts point E to larger 
T and smaller /i, and the vector-type interaction moves it in 
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Fig. 7: Phase (j> of the Polyakov loop as a function of (a) 6 and (b) T. 
Lattice data | 4, 8] are plotted by symbols. Curves represent results 
of PNJL calculations with set A. In panel (b), five cases (red, green, 
blue, pink and light blue) from top to bottom represent results of 
6/{tv/3) = 0, 0.4, 0.8, 1.0 and 1.2, respectively. The pink line 
terminates at T = Trw ~ 1.09 Tc, since is singular at T > Tbw 
in the case of 6 = n /3. 



the opposite direction. On the red solid curve between point 
E and point D both the first-order chiral and deconfinement 
phase transitions take place simultaneously. The light-blue 
dot-dashed curve moving up from point I represents the RW 
phase transition of first order, and point I is the critical end- 
point. The green dashed curve between point H and point E 
means the crossover chiral phase transition and the blue solid 
curve between point I and point E does the crossover decon- 
finement phase transition. Point F (G) is a crossing point be- 
tween the dashed (solid) curve and the ^ = line. Positions 
of points D-I are summarized in Table |lll] In panel (c), the 
pink dotted curve represents the lower bound of the location 
Me /Te of the critical endpoint E that the LQCD analyses of 
Ref. lUsf] predict. The position of point E in the case of param- 
eter set C is consistent with the results of the LQCD analyses. 



IV. SUMMARY 

We have tested the reliability of the PNJL model, compar- 
ing the model result with lattice data in the imaginary chem- 
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set 


D 


E 


F 


G 


H 


I 


A 


(2.02, 0.00) 


(1.84, 0.72) 


(0.00, 1.25) 


(0.00, 1.00) 


(i7r/3 X 1.53, 1.53) 


(i7r/3 X 1.09, 1.09) 


B 


(1.68, 0.00) 


(1.02, 0.87) 


(0.00, 1.05) 


(0.00, 1.00) 


(i7r/3 X 1.24, 1.24) 


(i7r/3 X 1.09, 1.09) 


C 


(1.80, 0.00) 


(1.51,0.72) 


(0.00, 1.05) 


(0.00, 1.00) 


(i7r/3 X 1.13, 1.13) 


(z7r/3 X 1.07, 1.07) 



TABLE III: Positions of points D-I in ^-T plane. The positions of these points are normalized as {fi/Tc,T/Tc) with Tc = 173 MeV. 
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Fig. 8: Phase diagram in the real chemical potential region, (a), (b), 
and (c) are calculated with the parameter set A, B, and C, respec- 
tively. Cross symbols with error bars indicate the latticle data taken 
from Ref. Points D-I are explained in the text. 



ical potential (/ij = TO) region. In this test, the model pa- 
rameters except Gv are adjusted so as to reproduce the mea- 
sured pion mass and decay constant at T = /i = and lattice 
data ® mill] at T > and /i = 0. In this step the eight- 
quark interaction plays an important role to make Tc (a) closer 
to Tc{'P) as discussed in our previous work [34]. With the aid 
of this, the PNJL calculation with the eight-quark interaction 
but without the vector-type interaction well reproduces the lat- 
tice data iS d] at finite 6* on ^ and Tc{<l>), but not on Tda) 
particularly near 9 — tt/S fully. The strength of Gv is then 
fitted so as to reproduce the data on Tc{(j) near 9 = tt/3. The 
primary result of the lattice simulations is that Tc{<P) coin- 
cides with Tc(o'), within numerical eiTors, in the entire region 
of 61 The PNJL model with the eight-quark and vector- 
type interactions can reproduce this property. Therefore, we 
can expect that the PNJL model with this parameter set is re- 
liable also in the /iR region. 

Finally, we quantitatively predict the phase diagram in the 
/iR region by using the PNJL model with the parameter set 
mentioned above. The critical endpoint does not disappear in 
virtue of the eight-quark interaction, even if the vector-type 
interaction is taken into account. This is the primary result 
of the present work. The lattice calculations at nonzero /ii 
have small lattice size (8^ x 4) |0,[8fl. Therefore, it is highly 
expected that lattice simulations with larger size will be done 
in the /ij region. 
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